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Abstract 

Using the ‘quantum optical approach’ we propose a model of multiplicity distributions in 
high energy collisions based on squeezed coherent states. We show that the k-mode squeezed 
coherent state is the most general one in describing hadronic mulitiplicity distributions in 
particle collision processes, describing not only pp collisions but e + e~, op and diffractive 
collisions as well. The reason for this phenomenological fit has been gained by working out 
a microscopic theory in which the squeezed coherent sources arise naturally if one considers 
the Lorentz squeezing of hadrons and works in the covariant phase space formalism . 


1 INTRODUCTION 

Although Quantum Chromodynamics is widely believed to be the theory of Strong Interactions, 
very few experimental results support this claim. In particular the behaviour of QCD at small 
momentum transfer i.e low energies is not understood. This lack of understanding reflects itself in 
the fact that particle production in high energy collisions cannot be explained within QCD. Given 
the absence of a detailed dynamical theory of strong interactions , one can adopt a statistical 
outlook and try to forecast macroscopic behaviour of a strongly interacting system given only 
partial information about their internal states. Experimental information about hadronisation in 
high energy collisions comes from the observation of jets of hadrons and the distributions of the 
final state particles. By using analogies with quantum optical systems one can get information 
about the types of sources( Chaotic, coherent, etc) that are responsible for hadronic emission. 
Also, by using adapting another quantum optical effect such as the Hanbury-Brown Twiss effect 
one can study the size and lifetime of the emitting region. This information can then be used to 
put restrictions on the microscopic theory pursued from the quark-parton end [1]. 

The experimental quantities amenable to the quantum statistical approach are: the multiplicity 
Distribution of final state particles (PIONS)given by 

P n = — (1) 

G inel 

where a n n-pion cross-section, the number of particles produced per unit rapidity dN/dy , where 
v = MIS) is the rapidity which plays the role of time in pion counting experiments, the 
moments of P„ and the two pion correlations which are analogous to Hanbury Brown Twiss effect 
for pions in rapidity space. 
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In particular, the quantum optical models are based on the assumption that multiparticle 
production takes place in two stages. In the initial stage formation of an excited system (fireball) 
which consists of a number of well defined phase space cells or ’sources’ which then hadronize 
independently. In these models an ansatz is made about the statistical nature of these sources 
and the resulting multiplicity distributions are compared with data [2], [3]. Table 1. gives the 
comparison of various quantum optical models. 


Table 1: Comparison of Quantum Optical Models of Multiplicity Distributions 


Nature Of Source 

Density matrix 

Probability 

Two pion 

One Source 

(Coherent State Rep) 

Distribution 

Correlations 

P(a) = ^ exp 

fi n r 

Pnm ~~ (i+fi)n+i °nm 

Geometrical 

S 2 (0) = 2 

P(a) = 6 2 (a — a') 

|a'| 3n e _l °' 1 * 
Pnn ~ r(n+l) 

Poissonian 

g 2 ( o) = i 

P(a) = 

O - nl -e-KPAt-HO 

Pnn ~ (l + *n»+i e 

xL ( -Kff \ 

XLn Vn(Hn)/ 

Glauber- Lachs 

1 < S 2 (0) < 2. 

K sources 




Gaussiaui (Chaotic) 

n+fc-1)!/ n/k \ n 1 

Pnn — n l(fc+l)I M+n/fc'' (1+n/fc) 4 

Negative Binomial 


Coherent +Chaotic 

_ W *21 

pnn (l+pn//c) n + fc 

— rprf 

l + pn/* Tk-l( -7* \ 

Xel iL n \ ( 1+pn/fc) ' 

Perina-McGill 

» » » » 


2 The Phenomenological model 

Experimentally there exists a large class of data (up) and low mass diffractive data that have 
mulplicity distributions with sub-Poissonian Statistics. Thus we seek a more general distribution 
than the ones given in table 1. A clue as to the appropriate distribution is that charged pions 
occur in pairs Furthermore the most general Gaussian source characterised by Gaussian Wigner 
Function . These facts point to the use of Squeezed Coherent states. 

We find that the k-mode squeezed state |o,r >= |ai,ri > |c* 2 ,r 2 > • • • | ak,rk > 
characterised by the multiplicity distribution: 

Pn fc = n, fc P n< - ^n i = n (2) 


p * 
* n 


X 


e [-fcQ 2 (l+x)]( 1 _ x 2 )fc/2(£)n 

2 

M l m H n ^J(Vky)2 2m 

m!(n — 2m)! 


/q 2 (1+x) 3 \i 

y = ( 2x ' 

7 = V; 7m = (7 + 1) " • (7 + (ra - !)) ; 7o = 1 


4 



and the second order correlation function: 


9k( 0 ) = 1 + 


2 sm/i 4 (r) + (2a 2 + l)sm/t 2 (r) - smh(2r) 
k(a 2 + sinh?(r )) 2 


( 3 ) 


is the most general distribution that fits a wide range of data [4]. If r > 0 there are regions where 
m < 1 and the distribution is narrower than Poissonian. If r < 0; <7*(0) is always greater then 
1 showing distributions which are broader than Poissonian. 

Hadronic distributions in pp collisions show broader than Poisonnian multiplicity distributions 
with a long multiplicity tail, which gets broader and broader with the increase of energy. The 
k = 3 mode distribution for n = 13.6, x = -0.20 and n = 26.1, x = -0.35 respectively fit 
corresponding ISR (62.2 Gev) and UA5 (540 Gev) data, a for each of these is thus fixed. 

To fit neutrino induced collisions in which the distribution is super-Poissonian ((4p) < 1), k = 
3,x = 0.5 fit data well. e + e“ collisions are fit by the k=2 squeezed coherent distribution with r 
close to zero . (nearly Poissonian.) 


3 The Statistics confronts the Dynamics 

We would now like to conjecture on the reason for this success and find an overlap with dyamical 
models. We search for incoming states of the hadronic fireball which will give rise to SQUEEZED 
COHERENT DISTRIBUTION. The candidate dynamical model of hadrons, which we find is 
appropriate is the covariant phase space model for hadrons which is a revival of Feynman et. al’s 
relativistic harmonic oscillator model [5] Kim and Wigner pointed out that the covariant harmonic 
oscillator model is the natural language for a covariant description of phase-space [6], [7]. In this 
paper, we use the covariant phase space distribution description of relativistic extended particles 
to give a phenomenological description of multiplicity distributions in the high energy collisions 
of hadrons 

Wave functions without time-like oscillations can be constructed by using the unitary repre- 
sentations of the Poincare group and imposing a covariant condition[8]. In this model two quarks 
bound together by a relativistic harmonic oscillator potential mapped onto 0(3,1) invariant har- 
monic oscillator equation. The ground state wave function tpp in the Lorentz boosted (primed) 
frame is 

V>oV) = [e~ ir)K3 }Mx) ( 4 ) 

where K 3 boost generator along the z axis, 

k 3 = i( 4 + 4 > (5) 

and 7 ] = Tanh~ l (P) 

In ‘Quantum Optical’ language, using light- cone variables we have: 

u=(t + z)/V 2 v=(t-z)/y / 2 «') 

Then in the Lorentz-transformed frame: 

(j'u = eV Qv = e ~ V( lv 

u ■- e~ T, u v 1 = e v v (7) 
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Introducing creation and annihilation operators 


d d 

a u = u+ — — ci v = v + — 
ou ov 


aj, = u — 


d 


d 


t w 

K = v 


du v dv 

we find that the wave function Tp% = ip(u',v') is a two mode squeezed state. 

Tpo{u',v') = |0,/3 >= |0, 77 |0, —77 >✓ 

The excited state is given by: 

I n,0>= (a|*') n |0 , 77 > (a|,/) n |0, — r/ > 

The condition for absence of time-like oscillations in the hadronic rest-frame 

K - o! v )\n,/3 >= 0. 

The physical wave functions Eire 

tpn(u', v') = \n,P>= m ) l n “ m > V >«' l m > ~V >v 


m=0 


( 8 ) 


(9) 


( 10 ) 


(11) 


in the Fock-space representation 


(“>')= E E 

ni ,U2 m=0 


Where [9]: 


for n,m even and 


n 

m 


G ni , n _ m (7?)G„ J>m (-77)^“ 1 (u)V'° 2 (w) 


^ rn\n\ ^ ( tanh(r))^ 

— V V l / \ / V o / 


cosh(rj ) /v 2 

sinh 2 (Tj) J 


m+n 

2 


™ n l!il ( =* ) x 

' sinh^Tn) ' 


^ (2A)!(m/2 — A)!(n/2 — A)! 


G n , m = (-l)^- 3 / 2 ( 


mini . tanh(r}) 


[ n ~ 1 rn — 1 1 

t 2 ' 2 i 


x E 


cosh(r ) ) ^ 
( sln/Pfr}) ) 


)’ 


l-l 


(2A + l)!(m - 1/2 - A)!(n - 1/2 - A)! 


( 12 ) 


(13) 


(14) 


for n,m odd. G n , m is non zero for both n,m even or both n,m odd , thus excitations of quarks 
occur in pairs, and the Lorentz squeezed vacuum is a many particle state . The above suggests 
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the identification of Hadronic sources in terms of squeezed states. 
In the ‘fireball picture ’ the Wigner function of the source is , [10] 


W n (u,v,q u ,q v ) 


x 


^(e^u’+e ’ , v 2 +e’>g2) 

7 r 


£ f " ) (-1 )X-JeV - e~' , ql}L m [e~ T 'v 2 - e’,, 2 ] 

m = 0 \ 1 1 


The number distribution for 1 particles in the n th 


excited state. 


Pi = 


n 


E E 

li+l2=l m= 0 


n 

m 


pr n "to) p h*(-’i) 


( 15 ) 


(16) 


pm,sg 


X 


(m)!/ 2 ! ( tanh(-rj) )m _ h 
(cosh(—r))) 2l * +1 2 


ru , \ 2,(™- / 2)tT 

F{-T],l2,rn) cos ( 


) 


(17) 


where: 


minfa h)/2 

F{-7],l 2 ,m)= J2 

X—0 


( -4 NA 

V sin/i 2 (— 77) / 


(A)!(/ 2 -2A)!(m-Z 2 -2A)! 


(18) 


The cosine terms imply Pi vanishes when | m — 1 2 \ or \n — m — 1\ \ is odd. so that the excited each 
of oscillator modes is excited in pairs. 

If each pair is associated with a two quark bound state(pion), the excited state con- 
tains pair correlated pions!!! 


4 Results and Conclusion 

The picture emerging is as follows the distribution of the fireball results from the excitation of 
oscillator modes of the colliding hadrons. This excitation takes place in pairs. Modes de-excite 
statistically emitting 2 pairs of quarks which we identify as two pions The phase space distribution 
of the fireball: 

| < n,/?|n, —0 > | 2 = (27 r) J dudvW£(u,v,q u ,q v )W0,(u,v,q u ,q v ) (19) 

Probability of emission of m particles from two independent populations 1 and 2 corresponding 
to each of the incident hadrons, forming an overlapping distributions is given as: 

m 

E (20) 

77l'= 0 

Total probability distribution thus becomes a product of the probability distribution of four 
squeezed sources: 

P m = E PZW p Z(-n) p ™M)P»H(-v') 

T71 1 + — TU2 +7713 -f 7714 = m 
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For target Projectile collisions ft = 0 thus the probability of emitting n’ particles is: 


p« = EE 

p=0 m=0 


n 

P 


n 

m 


/ ■ n ; + n . . .. .. . \l\l/2 ^ .tQTlhtirf) ' n ; -|-n . . mfn £ 

x (-1) 2 (p!(n — m)\m\(n ~p)0 7 ~ T7Z\ ( n ) 2 ( _1 ) 2 


cosh{rj) 


* E 


fa 


smh%|)) A+ ^ 


ST (2p)Kp/ 2 - P)!(^ - p)K™/2 - A)!(^ - A)! 


As 0 increases the distribution gets broader . 

For Central Collisions (3 = ft and by plotting mP m vs. for different values of 0 we see that 
the distributions become wider and skew symmetric as the value of 0 becomes larger. This is 
consistent with the variation seen in experimental data. 

The total probability distribution for the two nucleon system for n pions is: 


k/2 k/2 

p"= e nw e no-tf 

£]«.=»» * J2 n < =n * 


( 21 ) 


where k=6 for nucleon-nucleon collisions , k=4 for tvk collisions and k=3 for wp collisions (with 77 
positive). 

We include final state interactions in a simple fashion by assuming that the effect of interaction 
is to add coherence into the final state. This is consistent with the fact that in particle collisions 
experimental data shows some amount of coherence, especially in the low energy region , among 
the emitted particles. With the resulting density matrix we obtain the mutiplicity distribution for 
a variety of collisions and compare to data. The distribution we get is: 

k/2 k/2 

P*= £ pi P^’ coherent (a > r]) £ n P%' coh£rent {a,-n) (22) 

^ni=n * * 

Where the average number of particles emitted by each mode is given by:nj = a 2 + Sinh?(r l /) 
Above distribution fits the CERN ISR 62.2 GeV and UA5 540 GeV data. The k=3 distribution is 
compared with up data. The data is well reproduced by the distribution. For e + e _ collisions we 
take k=2 because the intermedeate state is the virtual qq state formed by the colliding electron 
and positron. 


In terms of hadronic final states the LEP energy (fts = 100 GeV) is equivalent to the SPS 
energy (fts = 546 GeV ) as far as total mutiplicities are concerned, in so far as It e (LEP) 
« lf p ( SPS) «26 . 

For the same value of n much narrower distribution for e + e“ distributions than the pp distributions. 
This is consistent with recent LEP data [11]. 

We can make some predictions for higher energies such as those observed at the LHC and 
SSC. Since widening of the distributions is related to the squeezing parameter 77 the lorentz boost 
of the hadronic fireball, at C.M.S. energies of 20 TeV and above we have a large 0 value and 
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higher modes will be excited. The multiplicity distribution for ultra-high energies is very broad 
and skew-symmetric, plot nP n vs.= for pp collisions for n = 50 

We can also calculate the Bose-Einstein Correlations of pions in this model by using the two 
mode state. Ongoing work is in progress to establish the connection of this model with QCD 
using the light cone formalism [12]. In this formalism it is also easy to incorporate temperature 
dependence by using Thermal Squeezed Coherent states. These would be of interest in heavy ion 
collisions. 
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